Energy levels and lifetimes of Nd IV, Pm IV, Sm IV, and Eu IV 
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To address the shortage of experimental data for electron spectra of triply-ionized rare earth el- 
ements we have calculated energy levels and lifetimes of 4/ n+1 and 4f n 5d configurations of Nd IV 
(n=2), Pm IV (n=3), Sm IV (n=4), and Eu IV (n=5) using Hartree-Fock and configuration inter- 
action methods. To control the accuracy of our calculations we also performed similar calculations 
for Pr III, Nd III and Sm III, for which experimental data are available. The results are important, 
in particular, for physics of magnetic garnets. 

PACS numbers: PACS: 32.30.-r, 32.70.Cs, 31.15.Ne 
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I. INTRODUCTION 



Magnetic properties of garnets are mostly due to un- 
paired 4f-electrons of the rare earth elements which are 
present in the garnets as triply- ionized ions (see, e.g. Q). 
Our interest in magnetic garnets was initially driven by 
proposals to use gadolinium gallium garnet and gadolin- 
ium iron garnet to measure the electron electric dipolc 
moment (EDM) 0,0). In our earlier works we performed 
calculations to relate the electron EDM to measurable 
values It turned out, in particular, that the value 

of the EDM enhancement factor is sensitive to the en- 
ergy interval between the ground state configuration 4/ 7 
and excited configuration 4/ 6 5d of Gd 3+ . The lack of 
experimental data on the Gd 3+ spectra produced an ini- 
tially large uncertainty which we were able to overcome 
by performing ah initio calculations of the Gd 3+ energy 
levels @. 

While gadolinium garnets seem to be the best choice 
for the EDM measurements 0, 0, other garnets may 
also be useful for the measurements of the T,P-odd ef- 
fects (EDM, nuclear anapole moment, etc.). Therefore, 
knowing the spectra of triply-ionized rare earth elements 
would be useful for future analysis. Unfortunately, the 
corresponding experimental data is very poor. In the 
present paper we calculate energy levels of the two most 
important configurations 4/ n+1 and 4/™5d for all rare- 
earth atoms for which corresponding experimental data is 
not available. This includes Nd IV (n=2), Pm IV (n=3), 



'Electronic address: dzuba@newt.phys.unsw.edu.au 
t Electronic address: lusa^oiu)@lidTeTmT ~ 

t Electronic address: johnson@nd.edu URL: www.nd.edu/~johnson 



Sm IV (?i=4) and Eu IV (n=5), but excludes Gd IV for 
which calculations were done previously [||. 

In the present paper, as in the earlier one we con- 
sider free ions only. Effects of the garnet environment 
on the energy levels of rare earth ions were considered 
by Dorenbos Q and can probably be used to extrapolate 
the spectra obtained here to the garnet environment. In 
any case, the spectra of free rare earth ions have their 
own significance. 

We also perform calculations for doubly-ionized ions of 
rare earth elements Pr, Nd and Sm, which have the same 
electronic structure as triply ionized Nd, Pm and Eu, re- 
spectively. Since experimental spectra are available for 
the Pr III, Nd III and Sm III comparison of theoretical 
energies with experiment allowed us to control the accu- 
racy of our calculations. 

A comprehensive set of theoretical data for the 4/™ +1 
- 4/" 5c? transitions including comparisons with data rec- 
ommended by the National Institute of Standards and 
Technology (NIST) @ is given in Tables I-IX of the ac- 
companying EPAPS document 0. 



II. CALCULATION OF ENERGIES 

Calculations in this paper are similar to those of our 
previous paper We use our relativistic Hartree-Fock 
and configuration interaction (RCI) codes as well as a 
set of computer codes written by Cowan [Tol ] and freely 
available via the Internet . 

We start our calculations from the relativistic Hartree- 
Fock method. Since we are dealing with open-shell sys- 
tems with a fractionally occupied 4/ shell, the Hartree- 
Fock procedure needs to be further specified. We per- 
form calculations in the approximation, which means 
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that all atomic electrons contribute to the Hartree-Fock 
potential. The contribution of the fully-occupied 4/ 
shell is multiplied by a weighting factor q/2(2l + 1) 
were q is the actual occupation number and I = 3 is 
the angular momentum. When the 4/ orbital itself is 
calculated, the weighting factor is further reduced to 
(q — 1)/(2(2Z + 1) — 1). Although the calculations are 
relativistic and the 4/ 5 / 2 and 4/7/2 states are different, 
we use the same weighting factor for both states. We also 
need the 5d 3 / 2 and 5d 5 / 2 states for the 4f n ~ 1 5d configu- 
ration. These states are calculated in the field of a frozen 
core from which one 4/ electron is removed. This cor- 
responds to the weighting factor (q — l)/2(2i + 1) for 
the contribution of the 4/ electrons to the Hartree-Fock 
potential. 

We now have four basis states 4/ 5 / 2 , 4/7/2, 5d 3 / 2 and 
5d 5 / 2 for the configuration interaction (CI) calculations. 
A calculation with these basis states and the standard 
CI procedure underestimates the energy interval be- 
tween the 4/" and 4/ n_1 5<i configurations. As shown 
by Savukov et al. [12j, this underestimate is due to the 
core polarization by the field of external electrons. In 
Ref. > correlation corrections to the energy of an ex- 
ternal electron in Ce IV or Pr V in the 4/ or 5d state were 
calculated to third order in many-body perturbation the- 
ory. It was found that core polarization increased the 4/- 
5d energy interval by about 20000 cm -1 , bringing theo- 
retical energies into very good agreement with the exper- 
iment. Similar energy shift of about 18000 cm -1 between 
the 4/ 7 and 4/ 6 5c? configurations of Gd IV were intro- 
duced "by hand" into the calculations with the Cowan 
code @. 

We use a somewhat different approach in our RCI cal- 
culations, where we simulate core polarization by adding 
a polarization correction 



to the HF potential in which the 5d state is calculated. 
Here a is polarizability of an ion in the 4/ n configura- 
tion, a is a cut-off parameter introduced to remove the 
singularity in the origin. We use a = as and treat a 
as a fitting parameter. Its values for different ions were 
chosen to fit energy levels of doubly ionized ions and the 
same values of a were used for triply ionized ions with 
the same electronic structure. Fine tuning the energies 
leads to slightly different values of a for different ions. 
The value a = 0.5a^ was used when data for the fine 
tuning was insufficient. 

Core polarization modifies both the 4/ and 5d states 
of an external electron; therefore, it would be better to 
modify both. However, to obtain correct energy inter- 
vals it is sufficient to modify either 4/ or 5d. Since the 
4/ orbitals are used for both the self-consistent Hartree- 
Fock calculations and the CI calculations, it is simpler to 
modify the 5c? states only. As demonstrated previously, 
, this procedure works very well for Eu III and Gd IV. 
Therefore, we use it for other rare earth ions as well. 



There is one further fitting procedure in the RCI cal- 
culations. We scale Coulomb integrals of second multi- 
polarity by introducing a factor f 2 = 0.8 before all such 
integrals. This scaling imitates the effect of screening 
of the Coulomb interaction of valence electrons by core 
electrons The value of the factor was chosen to fit 
energy intervals within the 4/™ _1 5d configuration of dou- 
bly ionized ions, where experimental data was available. 
The same screening factor was then used for all other 
ions. 

Calculations using the Cowan atomic structure code 
1 1 1| , which is also a configuration-interaction method (al- 
though in its non-relativistic realization), are similar to 
the RCI calculations. We use parameters available in the 
Cowan code to fit the known energies. One of these pa- 
rameters f 2 scales the Coulomb integrals similar to what 
was done in the RCI calculations. We use the factor 
f 2 = 0.85, which is very close to the value (0.8) used 
in the RCI calculations. Another parameter fog = 0.9, 
which scales the L — S interaction, is chosen to fit the 
fine structure of doubly ionized ions. Finally, the aver- 
age energy of the ground state configuration was shifted 
by hand to fit the energy intervals between the 4/™ and 
the 4/™ _1 5<i configurations. Note that this extra energy 
shift is due to correlation between core and valence elec- 
trons. The value of this shift is different for doubly and 
triply ionized ions. It follows that the Cowan code alone 
cannot be used to make predictions of the energy inter- 
val between the 4/" and the 4f n ~ 1 5d configurations. On 
the other hand, the procedure based on formula works 
well with the same parameters for both ions. Actually, 
energy shifts for triply-ionized ions used in the Cowan 
code were chosen to fit the results of the RCI calcula- 
tions. These shifts are 7000 cm" 1 for Nd IV, 9000 cm" 1 
for Pm IV, 10000 cm -1 for Sm IV, and 12000 cm" 1 for 
Eu IV. 



III. RESULTS 

In Table we present energies and corresponding 
Lande g-factors for some states of the 4/™ +1 and Af n 5d 
configurations of Nd IV (n = 2), Pm IV (n = 3), Sm IV 
(n = 4), and Eu IV (n = 5). For the 4f n 5d states we also 
present lifetimes r calculated by the Cowan code. There 
are two sets of theoretical energies: RCI energies (E R ) 
and Cowan code energies (E c ). For the 4/" +1 states we 
also present experimental energies (E N ) from Ref. 0. 
Note, that there are no experimental data for the Af n 5d 
states of these ions. 

The Lande g-factors presented in the table are useful 
for correct identification of the states. We present g- 
factors as they are calculated by both RCI and the Cowan 
code, together with the non-relativistic values, i.e., those 
obtained assuming pure LS coupling. The identification 
of states is usually easy for low states; however, states 
higher in the spectrum are strongly mixed, making level 
identification difficult, if not impossible. 
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Lifetimes of the 4/™5<i states presented in the table are 
calculated with the Cowan code, considering all possible 
electric dipole transitions to lower states of the ground- 
state configuration 4/™ +1 . Note that these calculations 
do not include core polarization and correlation effects. 
However, as shown in our calculations for Eu III and 
Gd IV 6], these effects are very important. For example, 
they change the lifetimes of Eu III and Gd IV states 
calculated in the Cowan code by a factor of 3.6. Similar 
corrections should be expected for other rare-earth ions. 
Thus, we believe that the lifetimes of the 4/ n 5rf states of 
Nd IV, Pm IV, Sm IV, and Eu IV in the Table dare too 
small by a factor of 3 or 4. 

More detailed theoretical results, together with com- 
parisons with available experimental data, are presented 
in the accompanying EPAPS document 0. Tables I - IV 
of give more detailed data for Nd IV, Pm IV, Sm IV, 
and Eu IV, respectively. Tables V - VII of 9] give more 
detailed comparison with experiment for the ions Pr III, 
Nd III and Sm III. Experimental data for these ions is 
sufficient for a reliable control of the accuracy of the cal- 
culations. The doubly- ionized ion Pr III has the same 
electronic structure as triply-ionized Nd IV. Similarly, 
the electronic structure of the pairs Nd III, PM IV and 
of the pairs Sm III, Eu IV is identical. Our calculational 
scheme can be roughly described as choosing parameters 
for a doubly ionized ion and then running the programs 
with the same parameters for the triply-ionized ion with 
the same electronic structure. Unfortunately, since there 
are no experimental data for the Pm III ion, we have no 
similar control of calculations for SM IV. Our confidence 
in the results for the Sm IV ion is based on the fact that 
the same procedure is used for all ions and very little ad- 



justment of parameters is needed to move from one ion 
to another. 



IV. CONCLUSION 

In summary, systematic Hartree-Fock and relativistic 
configuration interaction studies of energy levels and life- 
times of triply-ionized rare-earth elements in the 4/ n+1 
and 4f n 5d configurations have been presented. The re- 
sults of calculations using two different sets of computer 
codes for Nd IV, Pm IV, Sm IV, and Eu IV are in good 
agreement with one another and with experiment for the 
levels of the ground state configuration 4f n+1 . Our con- 
fidence in the accuracy of energy levels of the excited 
configuration 4f n 5d is based on similar calculations for 
Pr III, Nd III, and Sm III, for which sufficient experi- 
mental data is available. 
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TABLE I: Energies (cm x ), Lande ^-factors, and lifetimes r (sec) calculated by RCI (E n ) and Cowen code (E c ) 
= 2), Pm IV (n = 3), Sm IV (n = 4), and Eu IV (n = 5). Comparison with recommended NIST data, _B N by 
is presented. 
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